Mutual Information and Boson Radius in c = 1 Critical Systems in One Dimension 



Shunsukc Furukawa.^ Vincent Pasquicr,^ and Jun'ichi Shiraislii'^ 

^Condensed Matter Theory Laboratory, RIKEN, Wako, Saitama 351-0198, Japan 
"^Institut de Physique Theorique, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France 
'^Graduate School of Mathematical Science, University of Tokyo, Komaha, Meguro-ku, Tokyo 153-8914, Japan 

(Dated: March 5, 2009) 

We study the generic scaling properties of the mutual information between two disjoint intervals, 
in a class of one-dimensional quantum critical systems described by the c = 1 bosonic field theory. A 
numerical analysis of a spin-chain model reveals that the mutual information is scale-invariant and 
depends directly on the boson radius. We interpret the results in terms of correlation functions of 
branch-point twist fields. The present study provides a new way to determine the boson radius, and 
furthermore demonstrates the power of the mutual information to extract more refined information 
of conformal field theory than the central charge. 
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Conformal field theory (CFT) provides a power- 
ful framework to study one-dimensional (ID) quantum 
many-body systems. One can recast the low-energy 
degrees of freedom of both the bosonic and fermionic 
gases into a simple bosonic field theory, known as the 
Tomonaga-Luttinger liquid (TLL) |j. Thanks to the 
recent technological advances, a precise correspondence 
between the TLL predictions and various (quasi-) ID sys- 
tems, such as carbon nanotubes antiferromagnetic 
chains and cold Bose gases in ID traps II IS cur- 
rently being investigated. 

Given a microscopic model, an important and often 
nontrivial issue is how to obtain the effective field the- 
ory controlling its long-distance behavior. The notion 
of quantum entanglement, or more specifically, the en- 
tanglement entropy, has been extensively applied as a 
new way to address this basic matter. From a quan- 
tum ground state j^*), one constructs the reduced density 
matrix pA Tr^ |\E')(4'| on a subsystem A by trac- 
ing out the exterior A. The entanglement entropy is 
defined as Sa — Tr pA^ogpA- In ID quantum crit- 
ical systems, the entanglement entropy for an interval 
A = [xijXjJ embedded in a chain exhibits a universal 
scaling [10 I §,00: 

SA = ic/3)\ogix2-xi) + si, (1) 

where c is the central charge of the CFT and si is a non- 
universal constant related to the ultra-violet (UV) cutoff. 
This scaling allows to determine the universal number c 
as a representative of the ground state structure, with- 
out having to worry about the precise correspondence 
between the microscopic model and the field theory. 

As it is well known, the central charge is not the only 
important number specifying a CFT. In the bosonic field 
theory with c = 1, the boson compactification radius R 
(or equivalently, the TLL parameter K = l/(47ri?^)) is a 
dimensionless parameter which changes continuously in 
a phase and controls the power-law behavior of various 
physical quantities. It is natural to ask how to identify the 



boson radius as a generic structure of the ground state. 
In this Letter, we demonstrate that the entanglement 
entropy can achieve this task if we consider two disjoint 
intervals, A ~ [xi,X2] and B = [2:3, 2:4]. We analyze the 
scaling of the mutual information defined as 

Ia:B ■= Sa + Sb - Saub- (2) 

This measures the amount of information shared by two 
subsystems jlj, |l^ . A numerical analysis of a spin-chain 
model reveals a robust relation between Ia-.b and R, ir- 
respective of microscopic details. We compare the result 
with the general prediction of Calabrese and Cardy (CC) 
, and find a relevant correction to their result. 
Roughly speaking, the mutual information (H) may 
be regarded as a region-region correlator. It is known 
that Ia:B is non-negative, and becomes zero iff paub = 
PA ® Pb, i-e., in a situation of no correlation [0. A 
motivation to consider Ia-.b comes from that microscopic 
details at short-range scales, which are often obstacles 
when analyzing point-point correlators, can be smoothed 
out over regions. As we enlarge the region sizes, we ex- 
pect that Ia:B detects essential features of the correla- 
tions emerging in the coarse-grained limit. When there 
is a long-range order in local operators, we have Ia-.b 7^ 
for finite local regions A and B, even in the limit of large 
separation jl^. In a critical system with power-law de- 
caying correlations, Ia-.b goes to zero if A and B are far 
apart in comparison with their lengths, rA and rg. How- 
ever, if rA and rs are of the order of the separation, Ia-.b 
can remain finite, which is the situation we examine here. 

First, suppose we treat the mutual information (||) fol- 
lowing the prediction of Calabrese and Cardy Q , which 
will turn out in our analysis to correspond to the SU (2)- 
symmetric case. For an infinite chain, the entanglement 
entropy on double intervals AUB was predicted to be H 

c / a;2iX32a;43a;4i\ , , 

Saub = o^ogl + 2si, (3) 

with Xij = Xi — Xj. Here the constant term 2si is deter- 
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FIG. 1: (color online) The mutual information for fixed divi- 
sions rA:rc-rB:rD~l-l-l-l and 1:2:1:2, versus rj — 2tvR^ 



set the magnetization at M = with fc = 0, 1, 
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for -1 < A < 1 and with fc = 1, . . . , f - 3 for 1 < A, 
so that the system is inside the critical phase. Black and 
green points correspond to the larger (L = 28, 30) and smaller 
{L = 24) systems, respectively. Horizontal red lines indicate 
the Calabrese-Cardy result 



mined so that Saub ^ + Sb in the limit X2i,X4^ <C 
2^311 For a finite chain of length L, one replaces Xjj by 
the cord distance ^sin^ in Eqs. (|l|) and (|) g 
We now consider a division (r^, r^, td) of a finite 
chain as depicted in the inset of Fig. ^ Then the CC 
formula for the mutual information reads 



rCC 



3 log 



sin 



n{7-A+rc) „• T{rB+rc) 
L '""^ L 



sin ^ sin ^ 



(4) 



Note that the UV-divergent constant si has been can- 
celled out in the mutual information, and the resultant 
is invariant under global scale transformations. Sim- 
ilar ideas of eliminating the UV-divergence have been 
suggested by Casini and Huerta jl^ and have also been 
exploited in the context of topological entropy in 
higher dimensions. Henceforth, lengths of (sub)systems 
are measured in units of the lattice spacing. 

Now we turn to numerical analyses of the mutual infor- 
mation in a spin chain, based on Lanczos diagonalization 
of finite systems up to L = 30. We consider a spin-i 
XXZ chain in a magnetic field, 



H 



hY^S^. (5) 



3 5*1, is a 



Since the magnetization per site, M := j^"^ 
conserved quantity, we can label the ground-states using 
M. A c = 1 critical phase extends over a wide region 
in A > — 1 B. The boson radius R controls the nature 
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FIG. 2: (color online) Mutual information Ia-.b as a func- 
tion of ^ for divisions (rA, rc, r^, ru) = {r, ^ — r,r, ^ — r). 
We set h — 0, and symbols with different shapes correspond 
to different A = —0.8, —0.6,0, 1. Filled and empty symbols 
correspond to L = 28 and 24, respectively. 



of critical correlations. Indeed, the exponents for the 
leading algebraic decay of spin correlations, (S'JS'J,) and 
(S^S^,) - AP, are given by ?/ := 2ttR'^ and min(l/?7,2), 
respectively. Henceforth we frequently use 77 := 2t:R^ 
instead of R. For — 1 < A < 1 and = 0, we have 77 = 1 — 
(I/tt) arccos a. For A > 1, the system is in a gapped Neel 
phase at = and enters the critical phase at a critical 
field with 77 = 2. For general h ^ 0, rj can be determined 
by numerically solving the integral equations obtained 
from the Bethe ansatz |l9|, When increasing h, 
rj monotonically increases ( — 1 < A < 0) or decreases 
(0 < A) to i at the saturation. Summaries of the value 
of 77 in the M-A and h-A phase diagrams can be found 
in e.g. Refs. || and ^ 

We first evaluate Ia-.b for fixed divisions 
{rA,rc,rB,rD) - 1,1,1) and -1(1,2,1,2). Figure |l| 
shows a plot of Ia:B against 77 = 2ttR^ for various {M, A) 
in the critical phase. Remarkably, the data points al- 
most form a single curve for each type of division. The 
collapse of a two-dimensional M-A plane onto these two 
curves strongly indicates a direct dependence of Ia-.b on 
R. Agreement with the CC formula (jj) can be observed 
only around ry = 1 (5t/(2)-symmetric case). One can 
also observe that Ia-.b is symmetric under 77 — > I/77, 
which might reflect the duality in the effective theory. 

In Fig. ^ we plot Ia-.b as a function of j- for divisions 
{vA, rc,rB,rD) = {r, ^ ~ r,r, ^ — r), in comparison with 
CC formula (^). For each A, the results from L = 28 and 
24 obey a single curve, indicating the scale invariance of 
Ia-.b- The curve for A = 1 agrees well with the CC 
formula (Q). In other cases, the curves run above the CC 
formula. We can confirm that Ia-.b approaches zero in 
the limit ■£ — > 0, as expected for systems without long- 
range order. If we subtract the CC formula (see black 
circles in Fig. we find that the curves are symmetric 
under 4^ 



and have maxima at ^ 
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FIG. 3: (color online) The deviation of the "Renyi" mu- 
tual information /^"^ from the CC result for divi- 
sions (r, — r, r, -I — r). Different symbols correspond to 
71 = 1,2,3,4. 
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FIG. 4: (color online) Ia-.b—Ia/b versus the cross ratio x given 
in (p^. All the divisions {rA,rc,rB,rD) with 3 < ta < rB 
and 3 < rc < td are examined. For L = 28 (24), there are 
totally 305 (152) possibilities of such divisions. Black and 
Green symbols correspond to L = 28 and 24, respectively. 



As an extension of the von Neumann entropy, we also 
consider the Renyi entropy (or alpha entropy) defined as 



R' 



(") 



-1 
n - 1 



log(Tr pl) 



(6) 



The von Neumann entropy Sa can be reached in the limit 
n ^ 1. Following Calabrese and Cardy one can derive 
the following expression (originally found in Ref. ||]) for 
a single interval A = [xi,X2] in an infinite chain: 



R 



(n) 



1 



6n 



loga;2i 



(7) 



where s„ is again a UV-divergent constant. Likewise, 
within CC argument, the translation from von Neumann 
to Renyi can be done via replacements | "'g^c, Si — > 

(n) 

Sn- We define the "Renyi" mutual information as := 



R 



(n) 



R 



(») 



R 



(n) 

AUB- 



In Fig. ^, we plot the deviation of the "Renyi" mu- 
tual information 1^% from the CC prediction I^'^b^'^ ^ 
— ^^^c log [cos^ ^] . In contrast to the von Neumann 
case n ^ 1, we observe some oscillating dependence on 
j- for n > 1. Similar oscillations have also been reported 
for the single-interval entropy in Ref. |T| In Fig. || the 
oscillations in the n = 3 and 4 cases occur around the rel- 
atively smooth curves in the von Neumann case. From 
this, it is expected that .^^"b — ^a'b"'' consists of a smooth 
component, which depends little on n, and an oscillating 
component, which shrinks in the limit n — > 1 . The former 
should be controlled by the continuum description. 

Let us now discuss the origins of the deviation from 
the Calabrese-Cardy prediction and the dependence on 
the boson radius. We follow the formulation based on 
branch-point twist fields proposed by Cardy et al. ||2^ 
First, we represent the moment Tr p\ as the partition 
function on a n-sheeted Riemann surface TZn PI- Then, 



we relate it to a correlation function of twist fields T and 
T with conformal dimensions A„ = A„ = ^ ('^ ~ 
p2t . For double intervals AuB = [xi,X2]U [xa, X4] in an 
infinite chain, we can write it down as 



Tr pl^s CC {T{x,)f{x2)T{x;)f{x,)). 



(8) 



The SL{2, C) covariancc property requires this four-point 
function to have the following form: 



X2lX'i2Xi7iXii J 



2A„ 



X3,lXi2 



a:^2ia;32a;43a;4i 



2A„ 



Fn{x,x;ri), 



(9) 



with Xi = Xi. Here, Fn{x,x;r]) is a function of the 
cross ratios x := ^^^^"^ and x := ^2i£i3 normalized 



3:213:43 
3312:42 

as lim^^o Fn{x,x;r]) = 



3:213:43 
231842 ■ 

orn{x,x;rj) = 1, and should be determined 
by 77 = 2t:R^ as suggested by Fig The power func- 
tion part (. . . )^'^" (. . . )^^" in Eq. (|^ corresponds to the 
CC prediction and the function Fn gives an addi- 
tional contribution ^^^j log-F'„(x', x; 77) =: —fn{x;r}) to 

the Renyi entropy R'axjb- "^^^ mutual information de- 
tects this new part: 



;-(") 
1A:B 



^A:B 

tCC _ 
^A:B - 



= fnix;v), (10) 
lim fn{x;r]) =: f{x;r]). (11) 



The function fnix; 77) should satisfy (i) fnix; rj) ~> (x 
0), (ii) the crossing invariance fn{x;ri) = /„(1 — x;ri) 
required from R 



(n) 
AVJB 



^cuD f'^'^ ^ finite chain (see Eq. 
fn{x; l/rj) and/„(a;; 1) = 



(|1J) below), and (in) fn{x; 77) 
suggested by Fig. [|. 

As a check of this result, we plot Ia.b — I'a^b ^ ^ 

function of the cross ratio x in Fig. ^. For a finite chain, 
the cross ratio is given by 
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FIG. 5: (color online) n = 2 "Renyi" mutual information for 
the 1:1:1:1 division versus ri = 2-kR'^. The same symbols as 
in Fig. ^ are used. 

We can confirm that for a given A, and for various divi- 
sions (rAyfcTrByrD), the additional contribution to the 
CC result can be fit by a single curve with good accuracy, 
strongly supporting Eq. (pi]). 

For n — 2, two twist fields, T and T, are identical, 
and have conformal dimensions A2 = A2 = 1/16. The 
correlation function (^) of four twist fields with these 
dimensions (Ramond fields) has been derived previously 
p3|, p3. The crossing-invariant solution gives 



(13) 



where r is pure-imaginary, and is related to x via x = 
[()2{t)/63{t)]'^ . Here 62 and ^3 are Jacobi theta functions: 



^2 r 



(14) 



At a special point 77 = 1/2, Eq. ( p^ can be simpli- 
fied as f^{x;l/2) = log ^(1 + a;i/2)(i + (1 _ a;)i/2)/2. 
One can check that f^{x; rj) satisfies all the aforemen- 
tioned properties (i)-(iii). When two intervals of small 
lengths X21 — X43 =: r are separated far apart by a 
distance X31 — x^2 r), Eq. ( [T^ ) reduces to 

/I* ~ 2 (^^^^"""(^Z''''')^ This scales as the dominant spin 
correlation squared. 

In Fig. the formula f^ {x; rj) is drawn as smooth blue 
lines. The formula agrees relatively well with the data of 
I^^g — I^'^^'^^ for A = 0, and runs slightly above the data 
for A = -0.6. In Fig. |, we plot 1^^^^ for the 1:1:1:1 

division, in comparison with + I'^g^^ 
r(2) 



For 7? > 0.5, 

^^a-'b contains strong oscillations, and the formula goes 
inside these oscillations. For 77 < 0.5, oscillations are 
small, but the formula goes slightly above the data. This 
small disagreement might be due to finite-size effects, to a 
subtle difference between lattice systems and continuum 
descriptions, or to some missing factor in Eq. (|l3|). 



To summarize, we have shown that the scaling of the 
mutual information Ia-.b is controlled directly by the bo- 
son radius R. This result can be used as a new method 
to determine R from the ground state, complementary 
to the standard spectroscopic method based on the 
Drude weight and the compressibility. To obtain an ana- 
lytical expression of fn{x] rf) for general n and especially 
its n — > 1 limit /(x; rf) remains a challenging and intrigu- 
ing open problem. In general, we expect that every CFT 
has its characteristic function f{x) in the mutual infor- 
mation. This can be used as a fingerprint for distinguish- 
ing different CFTs, as originally suggested in Ref. 
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Note added. After the preprint of this paper was posted 
on arXiv, Calabrese and Cardy have added a note to their 
paper . Another analysis of the double- interval entropy 
has been done in parallel by Caraglio and Gliozzi ||2^ . 
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